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Abstract 

We discuss a toy model for an emergent non-relativistic gravitational 
theory. Within a certain class of Bose-Einstein condensates, it is possible 
to show that, in a suitable regime, a modified version of non-relativistic 
Newtonian gravity does effectively describes the low energy dynamics of 
the coupled system condensate/quasi-particles. 



1 Introduction 

In recent years the emergent gravity approach has constantly gained momen- 
tum pp. In this particular perspective, the gravitational field, encoded in 
the spacetime geometry and its dynamics, is seen as a kind of large num- 
ber/thermodynamical limit of some more fundamental theory. 

While there is no proof that this is the case, there are evidences to support 
this point of view. First of all, let us remind the striking correspondence existing 
between the laws of black hole mechanics in general relativity and the laws 
of thermodynamics [2]- This correspondence has inspired a deeper study of 
the thermodynamical aspects of gravitation. In particular, it has been shown 
that Einstein's equation can be seen as some sort of equation of state of a 
thermodynamical system at equilibrium [3|4]. Furthermore, from the study of 
the properties of gravitational collapse and the resulting black hole formation, 
very specific patterns have been recognized in the parameter space of the initial 
conditions which resemble the behavior of critical phenomena (phase transitions) 

mm- 



*To appear in the proceedings of the XXV Max Born Symposium, 28/6/2009-3/7/2009, 
Wroclaw, Poland 

t e-mail: sindoni@sissa.it 

'e-mail: girelli@physics.usyd.edu . au 

§e-mail: liberati9sissa.it 



1 



In addition, it has been realized that within several condensed matter sys- 
tems it is possible to distinguish some peculiar regimes in which the effective 
degrees of freedom are represented by fields propagating over effective pseudo- 
Riemannian structures [7]. 

For instance, in the case of perfect, irrotational and barotropic fluids, it 
can be proved that the perturbations in the velocity potential (i.e. the scalar 
function 9 whose gradient gives the velocity of the fluid, v oc V0) do obey a 
massless Klein-Gordon equation in a curved effective spacetime whose metric 
tensor is given by the so-called acoustic metric, 



_ P_ 



I -(c 2 s -v 2 ) : v l \ 



(1) 



where p is the local density of the fluid, c s it the (local) speed of sound and Vi 
is the velocity field of the fluid flow. 

Analogue models have been used to understand (and possibly to test in a 
laboratory) some peculiar aspects of physics in curved spacetimes, otherwise 
inaccessible (e.g. Hawking radiation). For a review of the subject see [7j. For 
the large majority, these analogue models for gravity do offer the possibility 
of studying some kinematical aspects of physics of curved spacetimes, leaving 
aside the issue of dynamics. 

Despite the fact that, to date, there are no completely successful models of 
emergent gravity, i.e. non-gravitational system possessing an effective General 
Relativity (GR) like limit, it is interesting to study simple toy models in which 
some aspects of the gravitational interaction can be mimicked. Here we present a 
toy model [8] , based on a well-known condensed matter system, a Bose-Einstein 
condensate of dilute, weakly interacting gas of bosons [§]• While this toy model 
is very far from representing an analogue of a realistic theory of gravity, it gives 
nonetheless some insights on some interesting features of emergent theories. 

The behaviour of dilute, weakly interacting Bose gases can be described 
with the formalism of second quantization, by introducing the field operators 
describing the gas atoms (V is the volume of the box in which we consider the 
gas) 

*( x ) = ^E a ^ k - x , (2) 

V k 

and a suitable Hamiltonian: 

H = |d 3 ^t (x) ^Z^!_ M+ ^ (x) |2^ (x) . (3) 

As it is well known, for p, k positive, this Hamiltonian leads to a condensation, 
i. e. the ground state of the system is characterized by a macroscopic occupation 
number of a single particle state. Within the formalism, this corresponds to the 
field operator '5 having a non-vanishing vacuum expectation value, 

$(x)«V(x)I + xto, (0|*(x)|fi)=V(x), (4) 
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where ip is a classical complex scalar field, the condensate wavefunction, de- 
scribing the mean field, and % is a quantum operator describing the residual 
quantum fluctuations around the condensate, or, in more physical terms, the 
atoms which are out of the condensate. 

Obviously, the nonlinearity in the Hamiltonian ^ makes the analysis of the 
physical properties of this system rather involved. As a first approximation, it 
is useful to assume that the non-condensed fraction is small compared to the 
condensed fraction. Consequently, one can introduce the parametrization 

* = ^ + £X + 0(e 2 ), (5) 

where e <C 1 is a dimensionless parameter which is introduced for bookkeeping, 
and organize an expansion in powers of e of the equation of motion. 



Gross Pitaevskii The lowest order equation describes the condensate wave- 
function alone, without the backreaction terms induced by the \ held, ft is 
called the Gross-Pitaevski equation and it is giving the first approximation 
to the condensate wavefunction dynamics. It is obtained from the equations 
of motion for the system by the replacement I n the setup we are 

considering, the equation reads: 

Bib h 2 

ifi^ = -^VV-M^ + ^lV, (6) 

which admits constant solutions, iji — (k/ij,) 1 / 2 , describing homogeneous con- 
densates. This kind of solutions, in the analogue model perspective, correspond 
to a flat acoustic metric, i.e. to Minkowski spacctime, and therefore will play 
a particular role in the following. The weak field limit of GR will correspond 
to the weak field limit around this constant configuration. In fact, as a direct 
calculation show, the correspondence between the acoustic metric and the 
complex scalar field ?/> is such that: 

, / fTi/^y + ,V • / \ => hoo oc u{x) (7) 
ip » (k/h) 1 ' 2 (1 +u(x) +iv(x)) J v ' w 

so that the perturbations in the number density encoded in u(x) are the natural 
candidate to represent some sort of Newtonian gravitational field. 



Quasiparticles The next order equation is an equation for the noncondensed 
fraction \ in the external field ip. The equation of motion has the shape: 

d h 2 

ih glX = ~^ v2 X + MX + MX t - (8) 

The Hamiltonian © does possess a global U(l) symmetry, ^> ~ ■» e Mx P, which 
is spontaneously broken by the condensation mechanism (leading to the ip ^ 
ground state). As a consequence of the Goldstone theorem, the excitations will 
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be gapless. (This can also be proved directly by diagonalizing the equation 
|J8]) using the Bogoliubov transformations.) Finally, it can be showed that, in 
some low energy regime, quasiparticles are described by a massless scalar field 
propagating on an effective acoustic geometry of the same shape of (Q]) . 

2 The coupled system 

In order to see how some sort of gravitational dynamics is encoded in the BEC, 
a suitable dynamical framework must be set up in order to see how the quasipar- 
ticles backreact over the condensate. This formalism consists in considering an 
improved version of the Gross-Pitaevski equation which consistently takes into 
account the effect of the particles out of the condensate. The Gross-Pitaevski 
(GP) equation is replaced by the so-called Bogoliubov-de Gennes (BdG) equa- 
tion 

d h 2 

ih—ip = -— V 2 -0 - urt) + K\ib\ 2 ib + 2nmb + nmip* , (9) 
at 2m 

where n, m are given by the expectation values: 

n = (E\x(x)^x(x)\E), m=(S\ x (x) 2 \E), (10) 

where the state |S) is the particular state one is considering. Notice that, if this 
state were the Fock vacuum state for particles, these expectation values would 
be identically zero. Notice also that one is implicitly taking a normal ordering 
in the particle operator, so that an unphysical (divergent) zero point energy is 
removed automatically. 

This equation, suitably modified, will give rise to the dynamics for the field 
u(x), and hence to an analogue of the Poisson equation. 

To construct some analogue of Newtonian gravity, we need massive particles 
as sources of the gravitational field (massless particles do not gravitate in New- 
tonian gravity). Therefore, the quasiparticles must not be Goldstone bosons, 
but instead pseudo-Goldstones: the U(l) symmetry has to be broken explicitly 
at the level of the Hamiltonian. In order to do this, we add an extra term to 
the Hamiltonian, 

H ^H = H + H\, H X = -±J d 3 x (*(x) 2 + (^(.t)) 2 ) , (11) 

where A is a coupling constant having the dimensions of an energy in these choice 
of units. The breaking of the U (1) symmetry has an obvious interpretation: the 
number of bosons is no longer a conserved charge. This could lead to some 
issues about the physical realizability of such a system. For a critical discussion 
on this point see [8]. 

Preliminary analysis of the condensate wavefunction based on the extension 
of the GP equation to this case, leads to the homogeneous solution 

i> = V^e- t9/ \ » C = ^, = 0, (12) 
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where the phase 6 = is fixed by stability of the condensate itself: for different 
values of 9 the quasiparticles would show a tachyonic instability. 

The analysis of the properties of the quasiparticles in the case of homoge- 
neous background (see [5] for details) leads to the conclusion that the quasipar- 
ticles dispersion relation is 

/ v 4 \ 1/2 

E =( Ma 4+M+h) ' (13) 

where 

In the case of small momenta, and when the condensate wavefunction is not ex- 
actly homogeneous (i.e. when u(x) ^ 0), the Hamiltonian for the quasiparticles 
takes the shape 



» j 2 ft 2 V 2 (m + A)(m + 4A) . . 
H g . P . « A^c 2 - W + ^ ^(x), (15) 



which leads to the identification of a "gravitational potential" : 

. (/i + 4A)(/i + 2A) 
*ffrw(x) = ^ £ w(x). (16) 

2.0.1 The emergent gravitational system 

Having presented the main ideas and required tools, we pass to the results. 
Consider the Hamiltonian with the U(l) breaking term. In the limit in which the 
backreaction of the condensate is small, i.e. in the limit in which there are few 
quasiparticles, when the condensate is almost homogeneous, the Bogoliubov-de 
Gennes equation can be rewritten as: 

V 2 --!) $ grav = 4nG N p ma ttcr + A, (17) 



where 



K (/i + 4A)(/i + 2A) 2 2k( ) u + 4A)(m + 2A) / 1 
GN = 47rtf m A3/^ + A)V2' A = W\ [ nn + 2 md 

Notice the peculiar splitting of the source term. A detailed analysis [8] shows 
that the expectation values (|10p always split into two contributions, one nonlocal 
term due to the quasiparticles (/? ma tter)j and an unavoidable vacuum contribu- 
tion, A, due to the inequivalence between the Fock vacuum for particles and the 
Fock vacuum for quasiparticles. 
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The reader will immediately realize that this interaction is very far away 
from a realistic Newtonian gravitational interaction: indeed, the would-be Pois- 
son equation includes a term which makes the interaction short ranged. In 
particular, this range is set by the healing length L, which is an UV scale for 
the physics of the BEC (very much like the Planck scale in quantum gravity). 
This might have been guessed from the beginning, since that the healing length 
represents the typical scale for the dynamics of the condensate. Despite the 
fact that this system is not an analogue for a realistic form of gravitational 
interaction, it does offer some intriguing inspirations which we develop now. 



3 Quasi-particles and locality 

An aspect of the physics of quasiparticles which is not often stressed concerns 
the issue of locality, i.e. whether the effective Lagrangian of the quasiparticles 
does obey the axioms of local quantum field theory. 

In BEC, the particles and quasiparticles field operators are related by a 
Bogoliubov transformation, which have the following general structure: 

uj A (k) = M£(k)cb B (k), (20) 

when working in momentum space. This relation is always linear. However, it 
is nonlocal: 

w A (x) = J(PyK£(x,y)<f> B (y), (21) 
where the kernel K is determined by the Bogoliubov coefficients: 

i^(x,y) = J d 3 kM£(k)e- lk < x -*K (22) 

This is the mathematical statement of the fact that quasi-particles are collective 
degrees of freedom. The structure of the transformation immediately implies 
that there is a mismatch between the notion of locality of the quasi-particlc with 
respect to the notion of locality of the atoms. As one easily realizes, the two 
classes of operators (f>, uj 7 separately, do obey canonical equal time commutation 
relations: 

[^(x),(0 B )t(y)] = ^ B J 3 (x-y), [^(x)> B )t(y)]=^ 3 (x-y), (23) 

which are a direct consequence of the fact that Bogoliubov transformations 
are preserving the algebra of the creation-annihilation operators. Therefore, as 
long as we use only one of the two families, there is no way in which a deviation 
from standard local quantum field theory can be manifest. However, the mixed 
commutators are nontrivial. It is straightforward to see that: 

[^(x),(^)t(y)] = (^)*(x,y), (24) 

Of course, this fact becomes crucial when the effective Lagrangian describing the 
physics of the quasi-particles involves terms mixing particle and quasi-particle 
operators. 
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(a) 



(b) 



Figure 1: (a) Feynman diagram of a four-particles vertex, (b) Lowest order 
Feynman diagram for the interaction of four quasi-particles (dashed lines) , in- 
duced by the interactions of particles (solid line). The circles denote the con- 
version of particles into quasi-particles at intermediate points. They represent 
the action of the Bogoliubov transformations. 

In the case of the BEC it is pretty clear how the underlying dynamics induces 
in the action for the quasi-particles an interaction term of the form 



which explicitly involves the particle fields, rather than the quasi-particles. 
Hence, when computing the effects of the interaction terms, e.g. scattering pro- 
cesses between the quasi-particles, the nonlocality encoded in the kernel K will 
necessarily enter the physical quantities (see figure [l}. 

The (effective) quantum field theories describing high energy physics are 
based on the assumption on locality, which is entering the physical observables 
at different levels. A detection of a violation of locality could be an indication 
that the (effective) quantum field theory we are using at low energies, namely 
the standard model of particle physics, is an effective field theory of quasi- 
particles over a condensate. Since we know that it is impossible to obtain the 
entire condensate dynamics out of a local quantum field theory based on quasi- 
particles (without knowing the Bogoliubov transformations), we could imagine 
that if the same picture holds for the real world, and in particular for gravitons, 
the difficulties in quantizing gravity are due to the fact that we are insisting 
using quasi-particle basis instead of using the particle one. 

It is interesting to note that the nonlocality is controlled by a very specific 
energy scale. We have seen that the Bogoliubov coefficients depend on the mo- 
mentum of the quasi-particles in such a way that quasi-particles operators are 
becoming particle operators for very high momenta. The scale of the crossover 
is given by the healing length scale, which represent then the typical scale as- 
sociated with nonlocality. 




(25) 
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4 Conclusions 



The model presented, obviously, is not a model of emergent gravity in the sense 
of GR. There are important conceptual and practical obstructions which forbid 
the BEC system to be used to discuss analogue of the gravitational field. In 
this sense, it is an attempt doomed to fail from the very beginning. First of 
all, the model is non-relativistic. Despite some encouraging results towards the 
extension to a relativistic theory made in [10 , it is still not clear how to evade 
the Weinberg- Witten theorem [Til H2 H3] in order to produce a viable model 
of emergent spin-2 gravity. Nevertheless, the brief discussion of the BEC-based 
toy model offers the possibility of showing very nice properties which could 
be of help in understanding some of the puzzles we find in the study of the 
gravitational field. 

The two main lessons which we should take concern the fate of locality 
and the smallness of the cosmological constant. As far as locality is concerned, 
there are two levels in a BEC: the particles together with their local interactions, 
and the quasiparticles which are the effective degrees of freedom propagating 
in the condensed phase. These latter are some collective degrees of freedom 
obtained from the non-local Bogoliubov transformations: the nonlocality of the 
relationship between particles and quasiparticles is translated into the failure of 
the interactions between quasiparticles to be described by a local quantum field 
theory. 

Furthermore, when the quasiparticles are seen as sources of the "gravita- 
tional field" encoded into the condensate wavefunction, the nonlocality mani- 
fests itself again in the non-local way in which they are coupled to it. 

A second important aspect which could be of interest in the discussion of 
the cosmological constant problem is the naturalness of the smallness of the 
cosmological constant term present in a BEC. As we have seen, the cosmological 
constant term is proportional to the depletion factor, i.e. by the ratio between 
the number of particles on excited states and the number of particles in the 
ground state (the condensate). Whenever this ratio is small, the cosmological 
constant is similarly suppressed. 

Notice that the smallness of the depletion factor depends on two properties: 
the dilution of the bosons and the smallness of the scattering lengths. 

The relation between the analogue of the cosmological constant and the 
depletion factor is intriguing. In fact, the smallness of the cosmological constant 
term is intrinsically related to the quality of the condensate phase, i.e. the 
larger is the condensation the smaller is the depletion factor. Accordingly, the 
smaller is the cosmological constant, the better is the mean field approximation, 
i. e. the way in which the condensate wavefunction encodes the dynamics of the 
condensate. Since the mean field is tied to the acoustic metric, the smallness 
of the "cosmological constant" term is directly related to how well the acoustic 
geometry will be able to describe the dynamics of the quasiparticles. 

This discussion should reinforce the idea that the emergent gravity program, 
and in particular the analysis of even simple toy models, could shed a light 
on very deep aspects of physics of spacetime beyond the grasp of standard 
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approaches like general relativity. 
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